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ABSTRACT 


The boundary condition for use with the diffusion solution for coupled radiation and 
conduction energy transfer is derived. An effective slip coefficient is presented as a 
function of the conduction-radiation parameter. Comparison to exact numerical solutions 
is good for the geometry of infinite parallel black plates. In the course of the analysis, a 
uniformly valid asymptotic expansion for the combined conduction and radiation problem 
in a nongray medium in the optically thick regime, and an exact analytical solution to the 
radiation-conduction transport equation in the boundary regime are obtained. 
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BOUNDARY CONDITIONS FOR THE DIFFUSION SOLUTION 
OF COUPLED CONDUCTION-RADIATION PROBLEMS 
by Marvin Goldstein and John R. Howell 
Lewis Research Center 

SUMMARY 

The diffusion solution for radiative energy transfer in gases contains the assumption 
that only local conditions contribute to the radiative flux at a point. When this assump- 
tion is justified, the diffusion solution is straightforward and accurate. For problems 
involving coupled radiation and conduction, the diffusion solution is probably the simplest 
approach that can be taken to determine temperature distributions. It is well known that 
diffusion methods for radiation are inaccurate near bounding surfaces, however, because 
energy from the boundaries, rather than the local gas volume, is a major factor. 

In this report, the correct boundary condition for the combined conduction- and 
radiation-diffusion solution is derived. The condition at the boundary is in the form of 
an effective temperature discontinuity that corrects for the error introduced in the dif- 
fusion solution by the presence of the boundary. This condition is a mathematical arti- 
fice only, because physically no discontinuity can exist in the presence of conduction. 

The effective slip is derived by the use of matched asymptotic expansions of the exact 
equation of transfer. The entire temperature distribution in the gas can be determined 
from the diffusion solution with a slip boundary condition and a correction near the wall 
found from the linearized exact equations. The linearized equations have been solved ex- 
actly by the methods of singular integral equations. Often the diffusion solution itself 
will provide the portion of the temperature distribution that is of interest. 


INTRODUCTION 

It is well known that the mathematical formulation of energy-transfer problems in 
enclosures containing absorbing- emitting gases in which radiation is the only mode con- 
sidered can lead to a prediction of a temperature discontinuity at the gas-boundary inter- 
face (refs. 1 to 3). The discontinuity arises because a gas element immediately adjacent 



to the boundary receives only a portion of the total absorbed energy from the boundary, 
while the remainder comes from surrounding gas and boundaries at different tempera- 
tures. The element then reaches an equilibrium temperature somewhere between the 
temperature of the adjacent boundary and the surrounding gas. This equilibrium temper- 
ature may be quite far from the boundary temperature. 

Heaslet and Warming (ref. 2) have presented predictions of the slip in the absence of 
conduction (radiative equilibrium) for a gray gas between infinite parallel plates. These 
authors used Chandrasekhar's tabulated X- and Y-functions (ref. 3) to provide useful ex- 
act solutions. 

Probstein (ref. 4) and Deissler (ref. 5) included the temperature discontinuity in 
formulating boundary conditions for use with diffusion solutions of radiative transfer. 

Although the temperature slip condition has proved useful in radiative transfer, it 
must be realized that an actual discontinuity in temperature arises only mathematically. 
Physically, some conduction of energy between the gas and the boundary will occur, ex- 
cept perhaps in certain extreme situations, and any conduction will remove the discon- 
tinuity. (In extreme cases, such as for an extremely rarified gas, other assumptions 
used in the equations of radiative transfer also become invalid, chiefly the assumption of 
local thermodynamic equilibrium. ) Thus, in any radiation problem involving a diffusion 
analysis or the typical transport analyses, no physical temperature discontinuity can ex- 
ist. However, for small conduction, very strong temperature gradients might be ex- 
pected near the boundary. 

The purpose of this report is to analyze the region of the gas near a bounding surface 
in which these strong gradients exist, and to determine the effect of thermal conduction 
on the gradients. An effective slip coefficient for use in the steady- state diffusion solu- 
tion for combined radiation and conduction problems is then derived for use with hot or 
cold boundaries. The utilization of the coefficients presented herein effectively extends 
the usefulness of the diffusion solution by specifying accurate boundary conditions for its 
formulation. 


SYMBOLS 


function arising from 

B v 

function given by eq. (87) 

separation of variable 

b 

intensity of emission in 

solution (eqs. (76) 


inner region 

and (77)) 


c v 

function defined by 

gas linear absorption 

V 

eq. (83) 

coefficient 


C-1 > Cn, Co, C A 

integration constants 

intensity of emission 
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D characteristic length 

E v function defined in eq. (103) 

G function defined by eq. (97) 

g function defined by eq. (98) 

H v function arising from separation of 
variables solution for h 
(eq. (76)) 

h >4 1) /T» 4 

I radiation intensity 

i intensity in inner region when sub- 
scripted 

K constant defined by eq. (31) 

k v function defined in eq. (103) 

L curve defined in fig. 2 

N radiation- conduction parameter, 

aX k /4aT^ 

q energy flux 

T absolute temperature 

t temperature in inner region 

u dummy variable of integration 

v separation variable 

X function defined by eqs. (101) 

x distance from bounding surface 

x + dimensionless coordinate, x*/6( e ) 

Z complex variable for eqs. (Ill) 

and those that follow 

a function defined by eq. (29) 

/3 defined by eq. (124a) 

r dimensionless energy flux, 

W T wl - T w2) 


y ratio of wall temperatures, 

"^w, 2 ^ T w, 1 

A temperature difference by extrap- 
olation of diffusion solution 
(fig- 1) 

6(e) intermediate or gage parameter 

e expansion parameter, 1/ aD 

0 function arising from separation of 

v (U 

variables solution for t^' 

(eq. (77)) 

9 angle measured from normal of 

boundary 

A k conductivity of gas 

li cos 9 

4 temperature gradient at x* — 0; 

dT*/dx*| x *. 0 

a Stefan- Boltzmann constant 

r optical thickness, ax 

$ function defined by eq. (93) 

^ radiation slip coefficient, 



£2 solid angle 

oo frequency 

Subscripts: 

c complementary solution 
D diffusion solution 

g gas 

k kinetic conduction 

R radiation 

r reference 

w surface or wall temperature 
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c<j frequency dependent 

1, 2 at surface 1 or 2 

+ > ~ 

approaching boundary from 
inside or outside domain 

Superscripts: 

* 

dimensionless property 

(0), (1), (2) terms in expansion 

— 

mean value 


ANALYSIS 

Consider a black wall at temperature T w that bounds a gas extending far from the 
wall (fig. 1). The mean free path for radiation 1/a^ is assumed to be much longer at 
all frequencies co than the gas kinetic mean free path, so that the gas can be assumed 
in local thermodynamic equilibrium. Scattering in the gas is neglected. 

The analysis, briefly outlined, is as follows: At regions far from the wall (i. e. , for 

O 

a^x — °°), a diffusion solution becomes exact to order (1/Da^) . At regions very near 



Distance norma! to boundary, x* 


Figure 1. - Definition of regions in gas near boundary. 
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the wall (a^x ~ 1), a linearized solution for the gas is shown to hold. The solutions for 
these two regions are then matched for intermediate distances from the wall. Up to 
terms of order (1/Da^)^ this expansion agrees with the usual diffusion solutions in re- 
gions not too near solid boundaries. On the basis of this correct expansion, the boundary 
conditions for the solutions to the diffusion equation must be modified by introducing an 
apparent discontinuity in temperature A at x = 0 to account for the boundary region 
where the diffusion solution breaks down. The values of A depend on the usual 
conduction-radiation parameter N and on the limiting value of the outer solution temper- 
ature gradient. 

The problem is now examined in detail. To determine the slip coefficient to first 
order, the only changes necessary to consider are those in the direction x normal to 
the wall. The results will be valid, however, for two- and three-dimensional problems. 

At any arbitrary position in the medium, the equation of transfer for radiation 
(ref. 3) can be written (scattering neglected) as 


cos 8 9 *oo 


Ct> 


ax 


= v t > - 


(la) 


where I is the local spectral radiation intensity, and B is the blackbody spectral 

UJ (jJ 

intensity. The equation denoting the conservation of energy between the conduction and 
radiation processes is, for conditions of steady state and no flow, 


d . dT 

A k 

dx dx 


_d_ 

dx 



I- , dco dfi 


At a solid black boundary, 



0 cos 9 < 1 


(lb) 


( 2 ) 


These equations can be nondimensionalized with the introduction of 
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/i = cos 6 





I* = 

CO 


CO 


4ctT: 


B* = 
co 


CO 


4crT' 


aX, 


N =• 


4aT! 


where a, D, and T r are arbitrarily chosen reference values of the absorption coeffi- 
cient, a characteristic length, and a temperature, respectively. These reference values 
are chosen in such a way as to make the nondimensional terms of order one (at least in 
regions far from the wall). Typically, a would be an appropriate frequency- averaged 
absorption coefficient, T f an average radiating temperature, and D a length of the 
same order as that over which the temperature changes significantly in the outer region. 
The general equations governing the conduction-radiation process then become 


M 

a* Da 

co 


3x* 


= B* 
co 


I* 

CO 


(3) 


and 


6 


I 


( 4 ) 


J_ d N dT* _ d 
aD dx* dx* dx* 


r f 

*oo=0 ‘ >/ is2=477 


til* d£2 dco 

CO 


Expansion for Outer Region 

For regions far from the wall, the analysis proceeds as follows: Let the quantity e 
be defined as 

_ _1_ 
aD 

By definition, an optically thick gas is one in which e is small as compared with one. 
Now equations (3) and (4) can be written 


01 * 

e _M ^ _ g* _ j* 


a* 9x* 
co 


CO CO 


( 5 ) 


dx* 


N f f p I* dco dJ2 

dx* dx* *oo= 0 *^2=477 


( 6 ) 


For e « 1, T*, I* , and B* can be expanded: 


I* = I (0) + elf 1} + e 2 I (2) 


CO (0 CO 


CO 


( 7 ) 


T*(x*) = T^°\x*) + eT^(x*) + e 2 T^(x*) + . . 


B*(T*) = B*? + eB f ( P + e 2 B, (2) + . 


co 


CO 


CO 


( 8 ) 

0) 


The relation for B* (T*) can also be expanded in a Taylor series about T 


( 0 ). 


B* (T*) = B* (T (0) ) + 


0B*(T^) 
co ' 

(oT" 


(T* _ t^^) + — 
2 ! 


0T 


j 3 2 B*(t(°)) 


(T* - T^) 2 + . . 


0(T^°V 


( 10 ) 
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Then, substituting equation (8) into equation (10) and gathering terms yield 


B* ( T *) = B* (T<°>) 


+ e 


3B* (T^) 
co ' 

jo r 


T ( : ) + e 2 


3T 


8B £( t(0) ) t ( 2) , i 82 b S< t(0) > r (i), 2 

( 0 ) 2 ! ^2 


3T V 


9 ( T (°)) 


+ . 


Comparison of equations (9) and (11) provides the relations 


( 11 ) 


bL 0) = B* (t(°>) 


,(D _ 


3 B *(t(°)) 


3T 


( 0 ) 


— T 


( 1 ) 


(12a) 

(12b) 


b ( 2) = 3B ^,( t(0) ) t (2) + 1_ 32b £( t(0) ) (1) 2 
w (0) 2! /nl 2 


3T 


g( X (°)) 


(12c) 


Equations (7) to (9) and (12) can now be substituted into equation (5). This substitution 
and equating like powers of e result in 


t( 0) = B (o) 

Cd Cd 


id) = B (i) . jl _ j®? 

“ “ a* 3x* 

Cl) 


gj(l) 

j(2) = b ( 2) . _y_ 

“ " a* dx* 

Cd 


(13a) 

(13b) 

(13c) 


3i( n ~ b 

j(n) _ B (n) _ _M_ 

“ " “ a* ax* 


(13n) 
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Similarly, substituting equations (7) to (9) and (12) into equation (6) and equating like 
powers of e give (after use of eq. (13)) 



Note that the zeroth power terms make no contribution because, since is independ- 

ent of ju, / B df2 = 0. Now define a dimensionless radiative conductivity as 
* 42=477 " 


A*(T*) = 



p 2 df2 


Then equation (14) becomes 


d 

dx* 


N 


dT 


( 0 ) 


dx* 


+ A.*(T' ') 


(0)x dT 


( 0 )" 


dx* 


= 0 


d 

dx* 



= 0 


(16a) 

(16b) 
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The term A*(T*) can be expanded about as was done for B* (T*) in equations (9) to 

(12), and the result is the expansion 


As before, 


Hence, 


A*(T*) = A*(T^) + e — — [A*(T^)]T^^ + 0(e 2 ) 

d T (°) 


A*(T*) = A^ + eA 


( 1 ) 


2 , 

+ e A 


( 2 ) 


+ 


x (°) _ a*(T^) 
x (l) _ dA*(T (0) ) T (l) 

d T (°) 


(16c) 

(16d) 


Multiplying equation (16b) by e and adding the result to equation (16a) give 


— Jn— (T^ + eT^) + A*(T^) — (T^ + eT^) + eT^ — 


dx* dx* 


dx* 


dx* 


A*( X (°)) 


= 0 (17a) 


Substituting the relation for A*(T*) into equation (17a) and using equation (16c) and (16d) 
result in 


dx* 


N — (T^ 
dx* 


+ eT 


(!)) + ( x(°) + eA^) — (t(°) + eT^) + 0(e 2 ) 
dx* 


= 0 


(17b) 


or using equation (8) results in 


dx* 


dT* , dT* 
N + A*(T*) 


dx* 


dx* 


+ 0(e ) = 0 


(17c) 


which is the same expression that is found by the usual diffusion approximation. The 

diffusion result can be obtained by expanding I* but not T*, and then neglecting terms 
2 

of order e . Such a procedure, however, does not allow the matching of inner and outer 
expansions to all orders in e to obtain a uniformly valid expansion, as is now shown. 


10 


r 


2 

However, equation (17c) is shown to be correct to order e in regions not too close to 
solid boundaries. Equations (13) and (16) are the governing equations for the outer re- 
gion. 


Expansion for the Inner Region 

Because the small parameter e multiplies the highest order derivative appearing 
in equations (5) and (6), the analysis deals with a singular perturbation problem wherein 
the formal asymptotic expressions (7) to (9) cannot satisfy all the boundary conditions of 
the original problem. This situation indicates the presence of a "boundary layer" in the 
region near the wall within which steep gradients are expected (i. e. , large values of 
91* /3x* and dT*/dx* in eqs. (5) and (6)). Thus the formal asymptotic expansion 
(eq. (17)) must break down in the boundary region, although it is expected to remain valid 
in the region far from the boundary. 

To obtain a valid expansion in the boundary region, the procedure continues with the 
introduction of a "stretched” coordinate (or optical depth) r: 

r = — = ax (18) 

e 


Introducing this new variable into equations (3) and (4) and maintaining the other dimen- 
sionless quantities the same as in the outer expansion give 


31 * 

JL w 

a* 3 t 

CO 


B* - I* 
co co 


(19) 


and 


A N ^ = i- f f n I* df2 dco (20) 

dT dT dT 0 «A2=4jr 

All terms in equations (19) and (20) are expected to be of order one in the inner region, 
since the highest order derivatives are no longer multiplied by e. When the coordinate 
is stretched, the gradients in the inner region are of order one. Therefore, an at- 
tempt can be made to expand the solutions to these equations in terms of the small pa- 
rameter e as 
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( 21 ) 


T*(t) = t^(r) + et^(r) + e^t^^(r) + . . . 


(22) 

(23) 

Expanding B* in a Taylor series about t^ and substituting equation (21), and then 
equating like powers of e between the resulting equation and equation (23) give 


t* -(0) -(I) 2.(2) 

I* = i' ' + £i v ' + e i v ' + . 
00 (0 0 ) 00 


B* = b<°) + cb W + e 2 bf) + 


b (0) = B * (t (0) } 
( 1 ) 


b 


3B*(t) 


o o 


at 


t(i> 


(24a) 


(24b) 


W°> 


b (2) = t( 2 ) 


oo 


9B*(t)“ 


a 2 B*(t) 

at 

t=t(°) 2! 

9(t) 2 


(24c) 


t=t«» 


Substituting equations (21) to (23) and (24) into equations (19) and (20) gives the linearized 
equations for the inner region as 


iL = - i< n ) 

a* dr w w 
00 


n - 0, 1, 2, 


(25) 


d M dt^ 
— N 

dT dr 


- f f 

dr ~'oo = 0 *A2=47t 



d£2 du> 


n = 0, 1, 2, . . . 


( 26 ) 


The boundary conditions (eq. (2)) at r = 0 must now be satisfied by the solutions to 
equations (25) and (26). By equating like powers of e in equations (21) and (22) at 
r = 0: 
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r 


t (0> 

ii 

H 

= 3 * 

(27a) 

t< n > = 0 

n > 0 

(27b) 

i(°) = B*(T ) 
CO O)' w' 

t-H 

VI 

=i 

VI 

o 

(27c) 

ii 

3 

V 

o 

0 < ii < 1 

(27d) 


Note that the conditions given by equations (27) are not sufficient to determine completely 
the solution of equations (2 5) and (26), and the conditions at large r must be further 
specified. 


Matching of Inner and Outer Solutions 


The formal expansions for the outer region (eqs. (7) to (9)) and those for the inner 
region (eqs. (21) to (23)) are still not completely determined. At this point, boundary 
conditions for the outer solutions cannot be specified because, as already noted, they 
break down near the boundary. The inner solutions are not specified completely because 
their behavior for large r has not been determined. This indeterminancy is removed 
by the requirement that the full expansion (inner and outer) be uniformly valid. In the 
present context, this uniformity means that, in a certain sense to be specified more pre- 
cisely, the inner and outer expansions must merge smoothly in some intermediate region, 
for example, 1/a « x « D. In this region, the procedure requires that both the inner 
expansions (eqs. (21) to (23)) and the outer expansions (eqs. (7) to (9)) be valid asymp- 
totic representations of the true solutions. 

For convenience, two order symbols are defined. First, 


<p(z) = 0[^(z)] 

means that there exists a positive such that 


limM 

z— 0 ['J>| 


<o9t 


and second, 
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<p ( Z ) = o[^(z)] 


means that 


limM = 0 
z— 0 |tf| 

The matching of the inner and outer expansions proceeds. It is well known (ref. 6) 
that to accomplish the matching a sufficient requirement is that an intermediate variable 


x 


+ 


x* 

6(e) 


exist with e « 6(e) « 1 and lim 6(e) = 0 (e. g. , 5(e) = e 3//4 ) such that, if j terms 

e— 0 

are retained in the expansions, the inner and outer expansions are expressed in terms of 
x + and all terms of o(e N ) are neglected with x + held fixed, then the two expansions be- 
come identical. That is, it is required that 


T( 0) [x + 6(e)] + . 


€ j T (i)[x + 6(e)] - t (0) 


lim^ 

£—0 


[x + e(£)l 


1 



- e J t 


L-(J) 




x + 6(e) 


-J V = 0 


(28) 


V. 


J 


with similar expressions for the expansions of I* and B* . The conditions on the inner 
and outer expansions must be found that will make them satisfy equation (28). 

First, define 


a*C T^ 0) ) s J A*(T^)dT^ 


(29) 


Then, substituting equation (29) into equation (16a) after integrating gives, for the outer 
region, 


Nt(°) + a*( T^) = Cj + C 2 x* 
From this equation, it is clear that 


(30) 
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( 31 ) 


lim T^[x + S(e)] = constant = K = T^(x* = 0) 
€-0 


From equations (12) and (13), it follows that 


lim I^ 0) [x + 6(e)] = B*(K) 


e— 0 


(32a) 


and 


l™ B™[* + e(£)] = B*(K) 
e— 0 


(32b) 


When noting equations (21) to (23) and that as e — 0, with x + fixed, t — it fol- 
lows that for the zero-order inner and outer expansions to match, 


lim t(°) 
e-0 


x + 5(e) 


lim t(°V) = K 


T’-^oo 


(33a) 


lim if 
e— 0 w 


x + 5(c) 


= lim 4 0 ) (t) = B*(K) 


T—°o 


(33b) 


lim b 
6-0 


( 0 ) 

CO 


x+ 5(e) 


= lim b^(r) = B* (K) 

7--*.oo 


(33c) 


To zero order, the inner equations (25) and (26) are 


9 i(°) 

JL = b (0) _ i(0) 

a* dr w w 

CO 


(34) 


and 


d_ 

dr 


N = -- f f jui^ dS2 dco 
dr dr -4>=0 - / S2=4 jt u 


(35) 


The zero-order inner equations (34) and (35), with boundary conditions (33) for r — °° 
and boundary conditions (27a) and (27c) for r = 0, have the solutions 
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(36a) 


K = T* 


+ ( 9 ) _ 

w 


b (°) = b* (T*) 
co co' w' 


(36b) 

(36c) 


b (0) = B* (T*) “I — M — 1 


CO CO' w 


(36d) 


-><0), . 


To match the first-order terms requires the expanding of a*(T' ') in a Taylor series 
around K as 


Q ,*| T ( 0 )[ x + 6(e)]j - o,*(k) + (T^ - K) 


da* 


dT 


( 0 ) 


+ O 


( T (°) _ K ) 2 


(37) 


T^=K 


Note that, at x* = 0, equation (30) becomes 

NK + a*(K) = Cj 


(38a) 


and also 

= A*(K) (38b) 

x (°) = k 

Substituting equation (37) into equation (30), using equations (38) and rearranging terms 
give 


da* 

dx (°) 


T 


( 0 ). 


C 9 6(e)x + 

K = — 

N + X*(K) 


+ O 


(T 


( 0 ) 



(39) 


Squaring equation (39) shows that 


(T^ - K) 2 = 0[6 2 (e)] 


(40) 


Now, 6(e) is chosen such that 
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(41) 


lim 
e— 0 


6 2 (e) 


= 0 


(For example, 6(e) = e^ where 1/2 < y < 1 satisfies this relation. ) In view of equa- 
tions (40) and (41), equation (39) becomes 


/q\ C 2 6(e)x + 

T (U) - K = — £ + o(e) 

N + X*(K) 


(42) 


After equation (30) is differentiated once and evaluated at x* = 0, 


C 2 = [N + A*(K)] 


dT 


(0) 


dx* 


x*=0 


(43) 


Substituting equation (43) into equation (42) gives 


- K = 6(e)x + 


dT 


( 0 ) 


dx* 


+ o(e) 


x*=0 


= 6(e)x H 


dT* 


dx* 


+ o(e) 


x*=0 


(44) 


Equation (16b) is integrated to obtain 


,(1) _ C 3 + C 4 X * 


N + X*(T^) N + X*(T^) 


When this is rewritten in terms of the intermediate variable and a new constant A* is 
defined by 


A* = T^(x* = 0) 


(45) 


the following equation is obtained: 
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T (1) [6(e)x + ] = 


N + x*(T^) 


+ 0[6(e)] 


= A* + 0[5(e)] 


( 46 ) 


The limit equation (28) is now written for j = 1 and equations (44), (46), and (36b) are 
substituted for , T^, and t^, respectively: 


r 


T (0) [x + 6(e)] + eT (1) [x + 6(e)] - t (0) 


lim< 
e— 0 


x + 8(e) 


- et 


( 1 ) 


x + 5(e) 


-A 


J 



K + 5(e)x 


+ dT 


(0) 


dx* 


+ o(e) 


x*=0 


+ e) A* + 0[6(e)]j - - et 


( 1 ) 


x*6 (e) 


= 0 


(47) 


Since K = T^, equation (47) becomes 


Urn 51 


( 0 ) 


e-0 e 


dx* 


- A* + t 


( 1 ) 


x*=0 


x+ 8(e) 


= lim<-^^ + 0[8(e)]l = 0 
e-0 e 1 


(48) 


Because r = x + 6(e)/e and 


lim(r) = lim 
e-0 e-0 

x ^constant 


x+ 5(e) 


(49) 


it can be concluded from equation (48) that 


lim t^(r) ~ A* + 
t — °° \dx*, 


x*=0 


( 50 ) 
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if the outer solution is to be correctly matched. Equation (50) then provides the boundary 
condition of the inner solution for t^ at large r. 

To determine the behavior of i^ at large r, expand equation (13a) in a Taylor 
series about T^ = T.* : 


4°)[x + 6(e)] = B<°>[x + 6(6)] = B^(K) + (T<°> - K)j 


9B^_ 


+ O 


( T (°) - K) 2 


(51) 


T^=K 


Then, the substitution of equations (40) and (44) into equation (51) gives (note that K = T.*) 


4 0) [ x+ 6(^)] = B^(T ; + ® 


(o r 


^ <0) ' T «>)_ T , 
w 


6(e)x 


+ dT* 


dx* 


+ o(e) 


x*=0 


(52) 


Substituting equation (12b) into equation (13b) gives 


3r(°) 3r(°) 

lLV»(0] = — *j- T (1 ) [x + 5(c)] - JL -f- 
(0) a* ox* 


9T' 


co 


(53) 


or noting that 


9B 


( 0 ) 


CO r + 


9T 


( 0 ) 


[x 5 (e)] = 


9B 


( 0 ) 


9T 


( 0 ) 


+ 0[6(e)l 


w 


HT 1 (®) + dT (®) 

51 — [x + 6(e)] = — — - 
dx* dx* 


+ 0[6(e)] 


x*=0 


dT* 

dx* 


+ 0[6(e)] 


x*=0 


and then substituting equation (46) into equation (53) to eliminate 


T (l) 


give 
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lW[x + 6(e)] = A*j 


3B^ 



3B^ 

00 


_ iL 

oo 

9T (°) 

'T'(®)_'T'^: 

W 

a* 

00 

dT (0) 


w 


dT* 

dx* 


+ 0[6(e)] 


(54) 


x*=0 


The matching condition for intensity is analogous to that for temperature (eq. (47)). Sub- 
stituting equations (54), (52), and (36d) into the matching condition gives 


lim< 
e— 0 


i 0) [x + 5(c)] + cl^x-tyc)] - I® L5M - eiW 5^] 


= lim^- 
e— 0 


1 ( 0 ) 


• 8(c)x + 2ll 

dx* 



9B® 

+ o(c) + € A* 


x*=0 

3T^ 


3B<°> 

€/i O) 


'T'-T* 


dT^ 

dx* 


x*=0 


" L « . 


r = 0 


(55) 


Note from equation (12a) that B* (T) = B^(T); thus, equation (55) becomes 


e— 0 


r 

• 







3B^ 


6(e)x + 

dT* 

4- A * - 

iL 

dT* 


9T (°) 

T - T * 

e 

dx* 

-f- 

x*=0 

a* 

00 

dx* 

x*=0 

V_ 

0 w 








id) 

oo 


x + 5(e) 


= limi^M + 0[6(e)] 

e— 0 e 


= 0 (56) 




or, finally, it can be concluded that 


lim i^V) 

J—CO 


A* + T - 


tJ- 


00 / 


dT* 


3B* 

00 

dx* 

x*=0 

3T* 


(57) 


T*=T* 

w 
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The boundary value problem for the first-order terms of the inner expansion is now com- 
pletely specified, and its parts are now collected in one place. From equations (25) and 
(26) with n = 1, equations (24b) and (27a) are 


Sid) 

j a to 

a* Sr 


St 


(t) 


h(l) 


id) 

to 


Jt=T* 


(58) 


and 


N dt^ = ^_ r r Mi (i) dS 2 do 

dr dr *4i>=0 ‘'£2=4tt 


_d_ Hf( J ) 
dr 


From equations (27b), (27d), (50), and (57), 


and 


i 1 ’ - o 


r = 0, 0 < q < 1 


td) = o 


T = 0 


tW-A. + r^) 


dx *4*=o 


X oo 


id) 

to 


A* + It - -tL 


dT* 


a* / dx* 

to/ 


x*=0 


SB* 

to 

ST* , 


l T*=T* 

W 


y — ► oo 


Now, set 


id)- id) + id) 

to to, p to, c 


(59) 


(60) 


(61) 


t ( MM» 


(62) 


where 
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I 



x*=0 



= A* + T 


dT* 

dx* 


id) = 

O), p 


A* + It - 


*w/ 


dT* 


3B* 

co 

dx* 

x*=0 

3T* 


W 


The direct substitution of expressions (61) and (62) into equations (58) to (60) then shows 
that id) and td) must satisfy the following boundary value problem: 


co, c 


aid) 

M co, c _ 

a* 3 r 

0) 


r 9B*\ 
, 3T* , 


hd) 


T*=T* 
W 


id) 

CO, c 


(63) 


dfd) /*°o 

N-^ = f f Mil?) 

dr *4o= 0 ~fo=4n- ’ 


dfi dco 


(64) 


For 0 < m — 1j 


i(!) _ f M dT* 

CO, c 


a* dx* 

CO 


- A*| 


'dB* 


co 


x*=0 


.dT*, 


T*=T* 
w 


t® = - a * 


T = 0 


T = 0 


(65) 


id) — 0 

CO, c 


*?> - ° 


T oo 


T — °o 


J 


( 66 ) 


Equation (63) can be integrated with respect to co to give 


M — 
dr 


r - Lid) do,= t ( i > f 

J as ,c c J 


f dB*\ 
co 

, 3T* , 


dco 


'Tijfc-'p* 

W 


£’c dw < 67 > 
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To proceed from this point, it is convenient, although not necessary, to assume that a* 
is independent of co; that is, the analysis deals with a gray gas. This approximation is 
made here only to eliminate a discussion of the effects of the spectral absorption coeffi- 
cient, which adds no new features. Under this assumption, a* = a* . Defining 


and using the expression 



equation (67) becomes 


( 68 ) 



a* dr 


rp 

T w ,(1) .(1) 

, c c 


(69) 


Equation (64), under the gray gas restriction, is 


dt 


(1) 


N 


dr 


= 2tt 


L 


1 i<»> 


I u dj u 


(70) 


and the boundary conditions (eqs. (65)) become 




JL dT * 

ia* dx* 


- A 


puts'} 

w 




0 ^ JLl < 1, T = 0 


x*=0 


(71) 


h(D- 


-A 


T = 0 


Defining 
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7T1 


l(l) 


h s. 


w 


t(l) 


T* 

w 


7* = a*T* 


I = • 


dT* 


T* a* dx* 
w 


x*=0 


. _ A* 
A = 

T 1 * 

W 


N 

w « 
w 


(72) 


equations (66), (69), (70), and (71) become 


0h f , 
p. = t - h 

(73) 

dT* 


N - 2 f hp dp 

w dr* S-l 

(74) 

h = pi; - A for 0 < jit < 1) 


It* - 0 

(75a) 

t = -A for 0 < ii < 1 J 


h — 0; t — 0 for t* — °° 

(75b) 


and solution of the boundary value problem specified by equations (73) to (75) will com- 
plete the solution of the entire problem. The method employed is similar to that used by 
Ferziger and Simmons (ref. 7) for a different radiation transport problem and is based 
on the work of Case (ref. 8). 
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Solution of the Boundary Value Problem 


A separation of variables solution is assumed valid for the problem specified by 
equations (73) to (75) and is of the form 


h = f H v 0a)e -T * /v A v dv (76) 

t= y A v 0 v e- r */v dv (77) 

Examination of the boundary condition, equation (75b), shows that 

R e (v) > 0 (78) 

The trivial solution h = t^ = constant exists, but it is neglected because of boundary 
condition (75b). 

Substitution of equations (76) and (77) into equations (73) and (74) gives 

"(v ■ ‘H = (79) 

N r l 

-y© v = 2 J H v (m)m dp (80) 

Eliminating 0 y from equation (79) by substitution in equation (80) and then rewriting 
equation (80) give 


(M - v)H y 


2v^ /*^ 

/ H (p)p dp 

N w J l 


( 81 ) 


and 


0 = - 


2v 


N 


w 



( 82 ) 
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In general, v may be a complex number. However, if v is assumed complex in 

equation (81), then for the imaginary part of v to be nonzero, N w must be negative. 

However, N is a parameter that can take on only positive values, so that for this prob- 
w 

lem Im(v) = 0, and v is taken as a positive real number by this argument and equa- 
tion (78). It is sufficient to restrict v to the range 0 < v < 1. With the function C y 
defined as 


C v' 




H y (p)p dp 


(83) 


equation (81) becomes 


2C v ! 

H y = ~ P- V. _J_ + 5(p - v)B y 


N 


(84) 


w 


li - v 


where P. V. denotes that the Cauchy principal values of the integrals are to be used. 
Multiplying equation (84) by p dp and integrating from -1 to 1 give 


c = 

V XT V 


N. 


W L 


2 + v In 


1 - v 
1 + v. 


+ vB 


(85) 


The function H v can be arbitrarily normalized so that 


N. 


C = v 
v 


w 


( 86 ) 


Equation (85) can then be written 


N. 


w r,. 2 3 1 - v 


B = _J2. - 2v - v u In 


1 + v. 


(87) 


With equation (87) B y is completely specified. 

Equations (84) and (82) become, with the substitution of equation (86)-, 


H = P- V. + B 6(p - v) 

P " v 


(88) 
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and 


0 


v 



(89) 


Now that H v and © v are known functions of v, they may be substituted into equa- 
tions (76) and (77) to yield 


h = P. V 


■/" 


3 -T*/v 


M - v 


A dv + B A e 
MM 


-t*/m 


"\ 

JLt > 0 


h = P. V. 


<4> 


> 


(90) 


dv 


M < 0 




and 


t - - f v 2 A dv 

4) v 


(91) 


The value of the remaining function A y is found by the substitution of equation (90) into 
the boundary condition (75a) to obtain 



This singular integral equation may be solved for A as follows: Consider the domain 
D (fig. 2) bounded by the closed curve L + L' in the complex plane. Define the function 
<f>(Z) as 
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y 



Figure 2. - Domain in complex plane bounded by closed 
curve. 


o 

As long as v°A is p-integrable (in the Lebesgue sense) on (0, 1) where p > 1, then 
<f>(Z) is analytic everywhere except on L. For such a function, the Plemelj formulas 
apply almost everywhere (refs. 9 to 11): 


and 


<$ + (v) - <f> - (v) = v 3 A v 


for 0 < v < 1 


<f> + (v) + <i> (v) = P. V. — 

u i 



u 3 A 


du 


u - v 


(94) 


(95) 


where e > 0 and 


$ ± (v) s lim $(v ± ie) 
e— 0 + 


Equations (94) and (95) are used in equation (92) to find that 

<f> + (v) = G(v)$”(v) + g(v) (96) 

where the function G(v) is defined as 
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G(v) s 


3 N - 4v 2 - 2v 3 
B v + jriv'* w 


In 


1 - v 
1 + v. 


7T1 


B. 


v 


7riv 3 .2 „ 3 

N w - 4v - 2v 


ln(- — ) + 7Ti 

.1 + V, 


(97) 


and 


t 4 3 A 

g(v) = - V - A - = 

B v - 7T iv 3 


2(gv^ - y 3 A) 


N w - 4v 2 - 2v 3 


ln( - — — | + iri 
1 + v, 


(98) 


Note that G(v) is continuous and nonvanishing on L and that g(v) is p-integrable in the 
range 0 < v < 1. The functions g^(v) and G^(v) on L + L’ are defined as 

g^ = 0, G^ = 1 on L’ 


gl = g, G i = G on L 

The behavior of G^ will now be examined. For all values of N , the modulus of 
G^(v) is 1. For the case of chief interest here, that is, for N w * 0, the argument of 
Gj(v) varies continuously in the range from 0 to less than 2u and back to 0 as v goes 
from 0 to 1 on L. The argument of G^ is 0 on L'; hence, Gj is single-valued on 
L + L ? , and the index of the boundary value problem (refs. 9 or 10) is 0. 

For the special case of N =0 (no conduction), however, the argument of G 1 (v) 
changes continuously from 0 to -2tt as v varies from 0 to 1; the index for this boundary 
value problem is -1. 

For an index of 0, a solution to the boundary value problem always exists for an ar- 
bitrary A and | (ref. 10). For an index of -1, however, an auxiliary condition must be 
imposed to relate A and This condition eliminates the possibility of satisfying the 
boundary condition T(x = 0) = Tq. The possibility of a physical slip in temperature at the 

wall for N =0 then arises. This case will not be pursued further. Note that 
w 


lim g(0 + e) = 0 
e— 0 

e > 0 

lim g(l - e) = 0 
e-0 

e > 0 
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and that and gj are therefore continuous and single valued on L + L f for + 0. 
The only discontinuity of <f>(Z) is on L so that equation (96) can be rewritten as 

<f* + (v) = Gj(v)$”(v) + g^(v) v contained in L + L’ (99) 

where the superscript (+) or (-) is understood to refer to the function approaching L + L’ 
from inside or outside D, respectively. 

The functions Gj(v) and g^(v) for any two arbitrary points except v=l satisfy the 
Holder condition 

|G(vj) - G(v 2 ) | < R[vj - v 2 |“ 

where R is a positive number and 0 < a < 1. Until recently, the Holder condition was 
required of both g(v) and G(v) before solution of equation (96) could be shown to exist. 
However, it has recently been established (refs. 10 to 12) that, as long as g(v) is 
p-integrable where p > 1 and G(v) is continuous as is the case herein, a solution for 
A in the class of p-integrable functions can be found that satisfies the boundary value 
problem almost everywhere (except on a set of measure zero); this solution is obtained 
(ref. 10) by setting 


where 


$ ± (v) = X ± (v) 



g(v) | 

X + (v) 


P. V. 



g(u)du 
X + (u)(u - v) 


X(Z) = exp 


1 

2ni 


f 1 ln °(Y) dv 

Jo v ' z 


( 100 ) 


(101a) 


and 


X + (v) - Gj (v)X - (v) 


X + (v) = exp 


- ln G(v) + P. V. 
2 


1 

277 i 



ln G(u)du 
u - v 






(101b) 
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With an appropriate choice of the branch of the logarithm, 


Define 



(102) 


(103) 


With the substitution of equation (100) into equation (94) and with the use of equation (101) 
and the definitions given by equations (102) and (103), the result may be written 



(104) 


If equation (91) is combined with equation (7 5a), the boundary condition becomes 


/ 


v^A v dv 


(105) 


Substituting equation (104) into equation (105) gives a relation between A/£ and N w> In- 
stead of proceeding directly in this manner, however, it is simpler to proceed as follows: 
Note that 


lim P.V. 

M— 0 



v°A /*1 o 

— dv = - / v A dv (106) 

H - v •'O v 
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2 

because lim v A is finite by equation (104). If the limit ju — 0 in equation (92) is 
v— 0 

taken, the following expression is obtained: 


-/ 


v A„ dv + lim BA = -A 


M 


M M 


(107) 


Comparison with the boundary condition (eq. (105)) gives 


lim BA =0 
p— 0 ^ ^ 


Substituting equation (104) now gives 


-A + lim E — P. V. 
v— 0 v 71 


du(i;u - A) 


/ 2 

E u( k u + 1 ) ( u - v ) 


= 0 


(108) 


Since (|u - A)/ 


uE u k u + 1 


is bounded at u = 0, and lim E is finite, the limit in equa- 

v-0 v 

tion (108) exists and reduces to 



du 


1 


S (k 2 

uy u 


uE„ k„ + 1 


1/2 E 


0 


/ 


du 


E, (k“ + 1 


u\ u 


1/2 


(109) 


where 


E n = lim E 
u v-0 


v 


Equation (109) provides a relation for A/£ as a function of N. This result can be 
further simplified as follows: 

Equation (101) is written as 


1 1 2?riv 3 1 

— = - — where veL 

X + (v) X~(v) B y - yriv 3 X + ( v ) 


( 110 ) 
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The coefficient of X + (v) on the right side goes to 0 at v = 0 and v = 1, and it follows 
from equation (101) that X(Z) never vanishes. 

The application of the Plemelj formulas and the boundary condition 


gives 



( 111 ) 


Expanding equation (111) in an asymptotic series near Z = °° yields 


X(Z) 


= 1 + 



dv 


E v k ; + 1 


1/2 


( 112 ) 


A similar expansion of equation (101a) gives 


- — = 1 + — — f In G(v)dv 

X(Z) Z 2rri y 0 


(113) 


and equating coefficients of 1/Z results in 


/ 

/ E 

n> v 


dv 


Since at v = 0, 1/ 


'( k v + ‘) 


/ 2 X 1 / 2 

1/2 


= — In G(v)dv = f tan — 
2i ^ 1 ^ 


— — )dv 

V 


E_ 


(114) 


is bounded, equation (111) becomes in the limit Z — 0 
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lim — — 
Z-0 X(Z) 



(115) 


It follows from equation (101a), since In G(v) = 0 at v = 0, that 


lim 


= lim exp 


Z-0 X(Z) u— 0 


— — P.V. f 
2ni J 0 


InGMdv 

v - u 


_ 1 _ 

E n 


(116) 


With the use of equations (114) to (116) in equation (109) finally results in 


A 




(117) 


From equations (8), (31), and (45), it follows that 


T*(x* = 0)-T^ = eA* + 0(e ) 


(118) 


that is, up to terms of 0(e ), e A*T r is the apparent discontinuity in temperature ob- 
tained by extrapolating the outer solution to the wall. 

As already shown by equation (17c), in the outer region the temperature satisfies the 

2 

usual diffusion approximation up to terms of 0(e ); hence, equation (118) gives the 
proper boundary condition for solutions to the diffusion equation, with the apparent dis- 
continuity in the wall temperature appearing to account for the region near the wall where 
the diffusion solution breaks down. If the radiation slip coefficient is defined as 






dv 


then equations (117), (118), and (72) show that 
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T*(x* = 0) - T* = eT* 4 U N w )| + 0( e 2 ) 

3 


= e 


1 dT* 


a* dx* 


- ^(N w ) + 0(e 2 ) 


x*=0 


or, written in terms of the dimensional quantities, 


3[T d (x = 0) - TJ 


4 1 


a dx 


WJ = <W w ) + 0(e 2 ) 


x=0 


where T n is the temperature obtained from the outer solution which is equivalent to the 

^ 2 
diffusion solution to order e . Because 


4T^ 3 [T*(x* = 0) - T*] = T* 4 (x* = 0) - T£ 4 + 0(e 2 ) 


w 


the expression for the slip coefficient can also be written 


4 4 

T£j(x = 0) - T 


w 


16 t 3 dT D 


3a w dx 


= i^(N w ) + 0(e ) 


x=0 


Alternatively, define as usual 


_ o dT^, 

■ 16 JL T 3 — 2 

3a w dx 


!x=0 


which is the radiation contribution to the wall heat flux in the diffusion approximation. 
Then 


T D< X = °) - < 




= ^(N w ) + 0(e z ) 


(119a) 
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APPLICATIONS 


Equation (117) provides a relation for the linearized radiation slip coefficient A/| 
in terms of a simple integral that depends on the single parameter N w , the usual 
conduction- radiation parameter. Examination of the equation (118) shows that the radia- 
tion slip coefficient ^ is related to A/£ by 


* = 



q r 


3 A 

4 £ 


(119b) 


where Tp(x = 0) is the temperature in the gas at the boundary to be used in the diffusion 
approximation. Deissler (ref. 5) predicts an approximate radiation slip on the basis of a 
second-order radiation diffusion solution, and his result for N w = 0 is 

^(N w = 0) = 0. 5 


The results from the present analysis, shown in figure 3, closely approach Deissler's 
result for small N . 


Prediction of Energy Transfer 

The slip coefficient derived herein can be readily incorporated into the diffusion solu- 
tion for predicting energy transfer. For the case of infinite parallel plates containing a 
gray gas and separated by a distance D, exact numerical solutions are available in the 
literature (ref. 13) for determining the accuracy of results. For this geometry, the dif- 
fusion solution (eq. (17c)) after one integration becomes 

q* = _n £2— - A*(T*) = — 2— (120) 

dx* dx* 4 ct t 4 

r 

in which the constant of integration is q* = q* + q* where q* is the dimensionless total 
energy flux and q* and q* are the dimensionless radiative and conductive fluxes, re- 
spectively. For no sources or sinks, q* is constant with x*, and equation (120) can be 
integrated directly. 

The slip coefficients ^ and are then used to eliminate the temperatures at 
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Slip coefficient. 



each boundary. From equations (119) and (120), 


( T *) 4 _ ( t *) 4 

*i 


= T- 


4q* 
1 + — N 


(121) 


w, 1 


where the sign of q* depends on the direction of the normal to the boundary and will be 
positive for i = 2 and negative for i = 1. Linearizing and substituting into the integrated 
equation (120) to eliminate TJ and T| give 


r = 




1 


— + +^n 

4 1 ^ 


1 + 


3N w .l* 1 -^ 

(i - r 4 ) _ 


(122) 


where y is the ratio of the wall temperatures T.* 2 / T w 1 and T D is the optical thick- 
ness of the gas aD. 

Predictions based on equation (122) are shown in table I and are compared with the 
numerical results of references 13 to 15 and with the simple additive solution of the exact 
uncoupled radiative and conductive fluxes as proposed by many authors. Agreement by 
both approximations is good, and the present formulation is slightly more accurate. An 


TABLE I. - COMPARISON OF APPROXIMATE AND EXACT SOLUTIONS FOR ENERGY 
TRANSFER THROUGH GRAY GAS BETWEEN INFINITE PARALLEL PLATES 


Optical 
thickness 
of diffusion 

Ratio of 
wall tem- 
peratures, 

Radiation- 

conduction 

parameter 

Exact 

Solution 

Present 

Additive 

solution, 

y 

at surface 1, 

Dimen- 

Refer- 

Dimen- 

Error, 

Dimen- 

Error, 

t D 


N 1 

sionless 

energy 

flux, 

r 

ence 

sionless 

energy 

flux, 

r 

percent 

sionless 

energy 

flux, 

r 

percent 

10 

0. 5 

1 

0. 336 

i 

0. 336 

0 

0.323 

-3.9 

10 

. 5 

. 1 

. 140 

i 

. 140 

0 

. 130 

-7.2 

10 

. 5 

.02916 

. 133 

14 

. 126 

-5. 3 

. 125 

-6. 0 

3 

. 2 

. 208 

. 583 

15 

. 567 

-2.9 

. 552 

-6. 5 

1 

. 5 

. 1 

. 863 

13 

.850 

-1. 5 

.773 

-10.4 

1 

. 5 

.01 

.647 

13 

.612 

-5.4 

. 581 

-10.2 
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II 


advantage of the present analysis is that temperature distributions in the medium are 
easily obtained by the integration of equation (120), but these cannot be obtained from the 
additive technique. 


More General Solutions to Optically Thick Situations 


For the purpose of discussion, consider the supposition that both the absorption coef- 
ficient a and the thermal conductivity are functions of the local gas temperature. 

A well-known fact is that, for a = a(T), pure radiation solutions of both energy transfer 
and temperature distributions may be greatly changed from the solutions for a assumed 
constant (ref. 16). For a = a(T), then, the additive solution for heat flux in combined 
radiation- conduction problems may be poor because of the influence of conduction on the 
temperature profile, which may change the radiative flux greatly from the radiative flux 

present with no conduction. 

o 

If terms of 0(e ) are neglected and the dimensional forms are substituted for all 
terms, equation (17c) can be written 


d_ 

dx 


MT) + 


16qT 3 

3a(T) 



= 0 


(123) 


Now define 


or 


a s 



A. (T)dT + — 

K o 



(124a) 




1 dB 
a(T) dT 


dT 


(124b) 


Equation (123) can then be written 


V 2 /3 = 0 (125) 

while the boundary condition (eq. (119)) is, with the use of the diffusion relation for q r , 
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(126) 


]3(x = 0) - j8 w 


i ** 

3 a ( T x=0 ) 9n 


x=0 


where d/3/dn is the rate of change in /3 with direction normal to the boundary. Equa- 
tions (125) and (126) are the general formulation of the diffusion equation for combined 
conduction and radiation with the extrapolated slip boundary condition. 


CONCLUDING REMARKS 

The correct boundary condition for use with the diffusion approximation for combined 
conduction- radiation problems was derived. The boundary condition is given in terms of 
a slip coefficient xjy, where 


* = 



<*r 


Values of i// as a function of the conduction-radiation parameter N are presented in 

w 

graphical form. Results obtained with the use of the diffusion approximation and the slip- 
coefficient boundary condition compare well with exact solutions for the case of a gray 
gas contained between infinite parallel black plates. 

The entire temperature distribution in the gas may be found by applying the solution 
of the linearized exact solution near the wall and the diffusion solution away from the wall. 
The procedure used herein guarantees that the solutions will match in the intermediate 
region and will give a uniformly valid representation of the entire temperature field. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, March 22, 1968, 

129-01-11-07-22. 
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